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Formation of VEGF isoform-specific spatial 
distributions governing angiogenesis: 
computational analysis 

Prakash Vempati 1 , Aleksander S Popel 1 and Feilim Mac Gabhann 2 * 
Abstract 

Background: The spatial distribution of vascular endothelial growth factor A (VEGF) is an important mediator of 
vascular patterning. Previous experimental studies in the mouse hindbrain and retina have suggested that VEGF 
alternative splicing, which controls the ability of VEGF to bind to heparan sulfate proteoglycans (HSPGs) in the 
extracellular matrix (ECM), plays a key role in controlling VEGF diffusion and gradients in tissues. Conversely, 
proteolysis notably by matrix metalloproteinases (MMPs), plays a critical role in pathological situations by releasing 
matrix-sequestered VEGF and modulating angiogenesis. However, computational models have predicted that HSPG 
binding alone does not affect VEGF localization or gradients at steady state. 

Results: Using a 3D molecular-detailed reaction-diffusion model of VEGF ligand-receptor kinetics and transport, we 
test alternate models of VEGF transport in the extracellular environment surrounding an endothelial sprout. We 
show that differences in localization between VEGF isoforms, as observed experimentally in the mouse hindbrain, 
as well as the ability of proteases to redistribute VEGF in pathological situations, are consistent with a model where 
VEGF is endogenously cleared or degraded in an isoform-specific manner. We use our predictions of the VEGF 
distribution to quantify a tip cell's receptor binding and gradient sensing capacity. A novel prediction is that 
neuropilin-1, despite functioning as a coreceptor to VEGF 165 -VEGFR2 binding, reduces the ability of a cell to gauge 
the relative steepness of the VEGF distribution. Comparing our model to available in vivo vascular patterning data 
suggests that vascular phenotypes are most consistently predicted at short range by the soluble fraction of the 
VEGF distributions, or at longer range by matrix-bound VEGF detected in a filopodia-dependent manner. 

Conclusions: Isoform-specific VEGF degradation provides a possible explanation for numerous examples of isoform 
specificity in VEGF patterning and examples of proteases relocation of VEGF upon release. 



Background 

Vascular endothelial growth factor A (henceforth called 
VEGF) is a critical pro-angiogenic factor secreted as 
numerous splice isoforms that together regulate the phe- 
notype and efficacy of growing vascular networks [1-7]. 
While the specific mechanism of this control is not fully 
understood, both isoform-specific receptor binding at the 
endothelial cell surface [8,9] and differences in the iso- 
forms' spatial patterning [6,7,10-13] are thought to be 
key. We have previously published studies of the impact 
of isoform-specific receptor binding [8,14], and here we 
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focus on the spatial patterning. The spatial distribution 
of VEGF isoforms is thought to be mediated by their 
interactions with heparan sulfate proteoglycans (HSPGs) 
in the extracellular matrix (ECM), and by proteases such 
as the matrix metalloproteinases (MMPs), which can 
cleave both VEGF [7] and the ECM [15]. Proteases have 
been shown to have important roles in inducing VEGF- 
mediated angiogenesis and tumorigenesis [7,16-20]. 
While it may seem intuitive that ECM binding regulates 
VEGF diffusion, computational studies suggest that at 
steady state, simple sequestration by HSPGs may have lit- 
tle effect on the soluble VEGF distribution. The specific 
mechanisms by which HSPG binding and proteolytic 
release regulate VEGF diffusion in vivo are therefore not 
yet fully understood, and we explore this here. 
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In mice, VEGF is primarily secreted as VEGF 12 o> 
VEGF 164 , and VEGF 188 (human VEGF is one amino acid 
longer: VEGF 12 i, VEGF 165 , VEGF 189 ) [3]; longer iso- 
forms include C-terminal motifs that increase binding to 
heparin and HSPGs in the ECM [21,22]. This increased 
matrix affinity can reduce the effective diffusivity of the 
isoform, altering the spatial gradient. For example, 
transgenic mice expressing only VEGF 12 o> an isoform 
lacking HSPG affinity [23], show a shallow VEGF gradi- 
ent in the developing hindbrain [13] and retina [12], 
whereas wildtype mice, which predominantly express 
the heparin-binding VEGF 164 , have a VEGF spatial dis- 
tribution that is markedly more localized (Figure 1A). 
Systems secreting VEGF 188 show the greatest levels of 
ECM and basement membrane VEGF deposition [6,24]. 

As noted above, VEGF sensed by endothelial VEGF 
receptors has two isoform-specific components: the spa- 
tial differences between VEGF isoforms; and the isoform 
specificity of binding to the receptors themselves as dis- 
cussed below. These give rise to a spectrum of isoform- 
dependent vascular phenotypes (Figure IB). VEGF regu- 
lates sprouting angiogenesis at the level of tip cell filopo- 
dial guidance and migration and stalk cell proliferation 
[12]. Systems (mice or tumor xenografts) secreting only 
the diffusible VEGF 12 o usually develop wide, tortuous, 
malformed vessels with infrequent branching (Figure IB) 
[6,12,13]. This phenotype suggests excessive proliferation 
and an insufficiency in sprout guidance. On the other 
hand, VEGF 188 alone gives rise to vessels that are thin 
and numerous and have high branching density [3,6]. 
VEGF 164 , similar to wild type mice, displays intermediate 
behaviour [4,6,13,25]. While these behaviours are 
thought to be primarily dependent upon the receptor tyr- 
osine kinases VEGFR2 (KDR/Flk-1) [12], VEGFR1 (Flt-1) 
also seems to play an important but not yet fully under- 
stood role, in part by modulating VEGFR2 signaling [26]. 



VEGF isoforms also differentially bind to neuropilin-1 
and -2 (NRP1 and NRP2), a class of semaphorin corecep- 
tors that function as coreceptors for VEGFR2. Neuropi- 
lin-1, for example, can greatly enhance the signaling of 
VEGF 165 relative to VEGF 121 [23,27]. 

Proteolytic release of VEGF from the matrix seems to 
biologically mimic secretion of shorter isoforms (Figure 
IB). For example, VEGF 120 elicits vascular effects at 
greater distances away from the source of secretion than 
VEGF 164 [6,12]; similarly, VEGF release by the protease 
MMP9 enhances VEGFR2 binding in quiescent vascula- 
ture by increasing soluble (diffusible) VEGF levels [16], 
activating tumorigenesis through the angiogenic switch 
[16,18,19]. Tumor xenografts expressing only the cleaved 
isoform VEGF 113 have large dilated vessels [7]. Inhibiting 
MMP9 results in markedly localized VEGF distributions 
reminiscent of VEGF 164 or VEGF 188 , with higher levels of 
matrix-bound VEGF [6,13,19]. Secretion of a synthetic 
VEGF 164 that is resistant to proteolysis (VEGF 164A108 _ 118 ) 
induces vascular patterning phenotypes similar to that of 
VEGF 188 alone, specifically high vascular densities [7] 
(Figure IB). While the direct cleavage of VEGF's C- term- 
inal heparin-binding and NRPl-binding domain is the 
most accepted mechanism of release (plasmin produces 
VEGF 110 from VEGF 165 [23]; MMP3 produces VEGF 113 
from VEGF 164 [7]), cleavage of HSPGs and ECM by pro- 
teases or heparinases is also a potent release mechanism 
[15,28-30], but one that liberates intact VEGF. Recently, 
degradation of soluble VEGF inhibitors such as connec- 
tive tissue growth factor and soluble VEGFR1 (sVEGFRl) 
has also been suggested to regulate VEGF [31-33]. 

The effects of protease-mediated VEGF release are not 
fully understood. VEGF 164A108 _ 118 - and VEGF 188 -expres- 
sing tumors show similar intratumoral vessel architectures; 
however, paradoxically, their ability to enhance tumor 
growth are vastly different: VEGF 188 either significantly 




Distance from secreting cells Protease- resistant VEGF164 Precleaved 

VEGF1 64 VEGF164 



Figure 1 VEGF isoform patterning and its effects on vascular patterning. VEGF is patterned in tissues in an isoform-dependent manner, 

with loss of heparin binding affinity translating to less localized, broader distributions [12,13]. A, Isoform specificity of VEGF distribution in the 

developing mouse hindbrain for mice transfected to express VEGF 120 versus wildtype mice, which express predominantly VEGF 164 [94]. B, 

Differences in VEGF localization are thought to give rise to vascular phenotypes that vary with the VEGF isoform, in a monotonic fashion. Larger 

isoforms, i.e. greater heparin binding affinity, give rise to sprouts that have greater filopodial directionality and vessels with greater branching 

density and smaller diameters. Similarly for VEGF proteolysis (red text): in tumors, loss of proteolysis exhibits similar features as VEGF 188 -secreting 

tumors, while VEGF 12 o-secreting tumors have similar features as tumors secreting precleaved VEGF [6,7]. 
I J 
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delays tumor growth [25] or cannot support it at all [6], 
while VEGF 164A108 _ 118 results in tumor hyperproliferation 
[7]. Similarly, while most studies implicate MMPs, e.g. 
MMP9, in pro-angiogenic and carcinogenic behaviors 
through VEGF [15,16,20,29,32,33], in others, VEGF pro- 
teolysis results in an angiogenic response that is uncoordi- 
nated and ineffective [17], leading to lower vessel density 
and decreased tumor growth [7]. 

It is not fully understood how the spatial patterning of 
different VEGF isoforms and the processing of VEGF by 
proteases control endothelial behavior, however several 
possibilities have been suggested by experiments 
[7,12,34]. For example, increased matrix-bound VEGF in 
tissues (due to expression of longer isoforms or inhibition 
of proteolysis) may increase branching behavior and 
decrease vessel diameters by increasing p38/MAPK sig- 
naling [7,35]. Alternatively, low levels of soluble VEGF 
may independently increase sprouting behavior [36]. On 
the other hand, vascular patterning might arise due to 
neuropilin-1 -dependent activation of p38/MAPK signal- 
ing [34], supported by the high affinity of VEGF 165 and 
VEGF 189 for NRP1 [37]. Finally, steep VEGF isoform gra- 
dients may improve tip cell filopodial stability while miti- 
gating stalk cell proliferation [12]. 

Whatever the dominant signaling modality that guides 
vascular patterning, it appears from experimental data to 
have two particular properties. First, the seemingly mono- 
tonic increase in vascular density and decrease in vessel 
diameter with VEGF isoform length (Figure IB) suggests 
that the critical feature of the VEGF distribution that con- 
trols vascular behavior also varies in a monotonic fashion. 
Monotonicity is also evident in the similar vascular pheno- 
types arising from mice dually expressing VEGF 120 and 
VEGF 188 , and those only expressing VEGF 164 , at the same 
total rate [6,13]. Second, increased VEGF cleavage has the 
opposite effect to increased HSPG binding [6,7] (Figure 
1A). All of the above modalities seem to support both con- 
ditions; for example, among the isoforms, VEGF 188 would 
be expected to have the highest levels of matrix-bound 
VEGF [24], lowest levels of soluble VEGF [21,22], most 
directional gradients [13], and possibly greatest NRP1 affi- 
nity [37]. Similarly, VEGF cleavage prevents VEGF binding 
to HSPGs and NRP1 [23], and by solubilizing VEGF 
should lead to greater diffusion. 

To separate out the above effects, we seek to develop a 
precise understanding of the effects of HSPG binding 
and VEGF proteolysis on the VEGF distribution by creat- 
ing a computational model. The model should replicate 
key experimental observations regarding the VEGF distri- 
bution (summarized in Table 1): longer isoforms with 
stronger heparin binding affinity have higher degrees of 
localization (shorter propagation distances and higher 
local peaks) (Figure 1A) [13]; increased heparin binding 
affinity results in greater levels of matrix deposition and 



Table 1 Goals of computational model to match with 



experimental observations 



VEGF Metric 


Effect of a decrease in 
isoform length 


Effect of VEGF 
cleavage 


Level of Soluble VEGF 


T [21] 


T [7] 


Level of Bound VEGF 


I [24] 


I [19] * 


Level of Total VEGF 


= [7] * 


= [7,19] * 


Magnitude of soluble 
VEGF gradients 


I [13] 


I [17] 


Spatial Range of soluble 
VEGF 


T [12] 


T [17] 



Note the concordance of the outcomes of decreased isoform length and VEGF 
cleavage. 

* hypothesized, or some supporting evidence, but no definitive experimental 
confirmation 



decreased levels of soluble-phase VEGF [7,21,22,24,38]; 
and VEGF-cleaving proteases should reverse these pat- 
terns [7,16,19]. The inverse relationship between soluble 
and matrix-bound VEGF gives rise to the question of 
whether total VEGF is conserved in vivo. While this is 
expected in vitro [22,39], and was also seen by certain 
studies in vivo [7,19], other studies found significantly 
higher total VEGF levels for VEGF 165 -expressing systems 
[25,39]. 

The current study extends our work on understanding 
the nature of VEGF transport in angiogenesis [40-44]. 
Among the numerous mechanisms of proteases modu- 
lating VEGF, we focus on VEGF cleavage, which we 
have previously studied [44]. Our model describes a 
sprout situated in a VEGF concentration gradient, to 
represent a range of biological systems, e.g. the hypoxic 
front of the hindbrain, tumor, or retina. 

We use the computational model (Figure 2) to under- 
stand mechanisms that might allow isoform-dependent 
VEGF localization and proteolytic redistribution to arise. 
Using a VEGF-VEGFR kinetic binding model, we investi- 
gate how VEGF signals may be detected by endothelial 
cells, e.g. soluble or matrix-bound VEGF, and the impact 
of NRP1. The results suggest that the system must exhi- 
bit isoform specificity in VEGF clearance or degradation: 
along with longer isoforms binding to the ECM with 
greater affinity, the longer isoforms must also be cleared 
or degraded from the soluble fraction at faster rates. 

Methods 

Model Formulation 

To study the role of different isoforms and VEGF pro- 
teolysis on the VEGF distribution, we constructed a 3D 
reaction-diffusion-based computational model of VEGF 
transport. The model is centered on a capillary sprout 
containing a single endothelial tip cell followed by stalk 
cells (Figure 2A), to study the VEGF concentration pro- 
file and gradients in the vicinity of the sprout, and to 
understand the receptor signaling and gradients detected 
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Figure 2 Schematic diagram of computational model. We simulate the transport of VEGF in an idealized microenvironment around a 
sprouting vessel. The above diagram, the base model, includes the processes of VEGF secretion and clearance (A), VEGF binding to sprout 
surface receptors (B), and VEGF binding to HSPGs and VEGF proteolysis in the interstitium (C). The domain is an axially-symmetric tissue cylinder, 
160 urn in length and 50 urn in radius, selected to mimic in vivo sprout spacing. The model considers the isoforms VEGF 12 i, VEGF 165 , and 
VEGF 189 in isolation to compare to single-isoform-expressing in vivo systems, and allows for cleavage of isoforms into VEGF 114 , the MMP-cleavage 
product. Receptor signaling is modeled with the isoforms' interactions with VEGFR1, VEGFR2, and NRP1 (B). While the reactions are shown 
primarily for VEGF 165 , isoforms VEGF 12 i and VEGF 189 , where bracketed, participate in similar reactions. VEGF and receptor-binding gradients are 
measured between the front and rear of the tip cell (indicated in gray arrowheads) (A). Extensions of this model explored here include gradients 
of HSPG and VEGF degradation (Figure 4). 



by VEGF receptors (primarily VEGFR2) on the endothe- 
lial tip cell. To dissect the behavior of individual iso- 
forms, we consider separate systems expressing a single 
isoform in isolation, similar to those tested experimen- 
tally in developmental and tumor models. 

A schematic of the VEGF transport and reaction is pre- 
sented in Figure 2B. In our model, we label our VEGF 
using the human isoforms nomenclature and represent 
cleaved VEGF by VEGF 114 , which is the predicted MMP- 
cleavage product of human VEGF [7]. While Figure 2 pri- 
marily details the interactions of secreted VEGF 165 , we 
also separately consider VEGF 12 i and VEGF 189 , and 
hypothetical isoforms of intermediate lengths and prop- 
erties. These intermediate isoforms are generated by 
smoothly varying the relevant parameters: HSPG affinity, 
NRP1 affinity, and ability to bind the VEGFR1-NRP1 
complex. While VEGF 121 has weak affinity to NRP1 [45], 
we assume it does not bind NRP1, and this does not 
change our results [43]. Intact (unproteolysed) VEGF 189 
may be unable to bind VEGFR2 [46], but binding of long 
isoforms of VEGF to VEGFR2 is predicted to be low in 
our model due to degradation and cleavage; cleaved 
VEGF represents the majority of VEGFR2 binding in this 
case. 

We assume that VEGF is secreted from the leading 
face of the cylindrical domain and is cleared by two 
mechanisms: internalization by the sprout's VEGF 



receptors; and clearance at the trailing edge of the 
domain (which physiologically represents the effects of 
both diffusion and internalization by downstream vascu- 
lature). As a base model (Figure 2, 3), we consider that 
the sprout expresses all three VEGF receptors, in a com- 
bination where excess NRP1 plays a potentiating role for 
VEGF 165 binding to VEGFR2 [27]. The sprout is situated 
in a cylindrical domain consisting of an acellular extra- 
cellular matrix (ECM) that contains uniform HSPG 
binding sites for VEGF. For VEGF proteolysis, we 
assume a uniform, single-step irreversible first-order 
MMP-mediated cleavage of VEGF to VEGF 114 (devel- 
oped previously in [44]). The dynamics of protease acti- 
vation and inhibition have been studied in the past 
[47-49], but are neglected here because of uncertainties 
in the transport of the protease itself. 

We have previously shown that the thin basement mem- 
brane layer does not significantly affect VEGF diffusion 
near the sprout [44], and thus do not consider the base- 
ment membrane here. As in previous models [41,42,50], 
receptors are considered pre-dimerized and binding of the 
ligand is a reversible single-step reaction; in addition, 
internalization and receptor insertion are balanced to keep 
plasma membrane receptor levels constant. Finally, though 
matrix-bound VEGF may be relevant to VEGFR2 activa- 
tion [35], our model only considers receptors binding to 
soluble VEGF. 
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Figure 3 Steady-state distribution of VEGF is not dependent upon interstitial HSPGs. Using the model, we calculated the steady-state 
distribution of VEGF 165 , in the presence (A) and absence (B) of interstitial HSPGs, assuming no interstitial proteases, and no other differences 
between the two systems. The VEGF distribution was initialized so as to have a mean soluble VEGF concentration of ~1 pM and a fractional 
gradient of ~5%/40 um at the leading edge of the sprout (z = 0); HSPG binding and receptor-mediated internalized were calculated. Conditions 
of the model are defined in Table 4: HSPG = 750 nM, VEGFR1 = VEGFR2 = 10 4 /cell, NRP1 = 3*10 4 /cell; secretion and clearance of the soluble 
fraction occurred at q = 5.27-1 0" 5 molec/um 2 -sec and k dear = 0.0399 um/sec. Traces of steepest descent arising from r = 0, 6, 20, and 38 urn and 
z = 80 um are overlaid to demonstrate the effect of internalization on the curvature of the VEGF distribution around the sprout. We also 
measured VEGFR1, VEGFR2, and NRP1 receptor occupancies, absolute and fractional gradients, and the sensitivity of relative gradient detection 
are given for each case. The sensitivity is defined as the ratio of the receptor gradient to the ligand gradient. 



We study several specific models with increasing com- 
plexity, to understand the mechanisms behind gradient 
formation. The model components as described above - 
diffusion, receptor binding, proteolytic processing - com- 
prise the "HSPG-binding-only" model This is so named to 
contrast with other models that add additional processes 
relevant to VEGF patterning in vivo: patterned HSPG gra- 
dients (the "HSPG-gradient" model) and VEGF degrada- 
tion. VEGF degradation is defined here as any interstitial 
process that destroys VEGF activity, without generating an 
active cleavage fragment as proteolysis does. As with pro- 
teolysis, degradation was characterized by a first-order rate 
constant, k deg . The study of degradation was done with 
two models: one assuming that degradation affects only 
soluble VEGF (the "soluble VEGF degradation" model); 
and one in which both soluble and matrix-bound VEGF 
are degraded (the "matrix-sequestered VEGF degradation" 
model). We will see later that the soluble VEGF degrada- 
tion model, in which HSPG-bound VEGF is protected 
from degradation, is effectively isoform-independent 
degradation; while degradation of HSPG-bound VEGF 
results in isoform-dependent degradation. Longer isoforms 
have higher degradation (despite having the same rate 
constant) due to longer retention times. 



In studying receptor binding, we incorporate a sprout 
expressing VEGFR1, VEGFR2, and NRP1. This sprout 
can itself influence VEGF transport (by receptor- 
mediated uptake), and to study the patterning of VEGF 
that would arise only by diffusion, matrix binding, pro- 
teolysis and degradation, we also considered a variant of 
the model that contains no sprout. 

VEGF/Protease Transport and Reactions 

VEGF biochemical reactions and transport processes are 
described by partial differential equations with appropriate 
boundary conditions. The computational domain consists 
of two components: the interstitium, into which VEGF is 
secreted, and a sprout surface layer, on which VEGF cap- 
ture by the receptors is modeled. Upon secretion from the 
leading edge of the domain (z = +L), VEGF transport 
within the domain is governed by the mass balance 
(VEGF 165 equation shown; similar equations for other iso- 
forms): 

= D 165 V 2 [V 165 ] - k^[V 165 ] - k v deg [V l65 ] 

k v,H ' (1) 

-T^[Vi 65 ][H] + k^ H [V 165 H] 
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Here [V 165 ] describes the volumetric "bulk" concentra- 
tion of the soluble fraction of VEGF 165 , assuming an iso- 
tropic porous matrix with available volume fraction 
K EC m- Di65 is the effective diffusion coefficient for 
VEGF 165 in the matrix; [H] is the density of VEGF-bind- 
ing sites, with reaction rate constants, k on v ' H and k 0 ff V ' H - 
k P describes first-order net rate of proteolysis, and k deg 
describes first-order degradation of soluble VEGF. We 
assume that k P and k deg are spatially uniform. The 
values of these parameters are given in Table 2 and 
described further below. 

Cleaved VEGF (VEGF 114 ) is generated from the proteo- 
lysis of free VEGF or matrix-bound VEGF (Figure 2C). 



3[Vii4 

at 



= D 114 V 2 [V 114 ]-fe d / eg [Vii4] 
+ k^[V 165 ] + k^ H [V 165 H] 



(2) 



VEGF proteolysis also produces an HSPG-binding C- 
terminal fragment [23] that maintains heparin affinity, 
but we do not include it explicitly as HSPGs are not 
saturated (i.e. [V]<< I<d V ' H ). The matrix components, 
HSPG and VEGF-HSPG do not diffuse: 



3[H] 

at 



l,V,H 



Kecm 



[V 16 5][H] + k^ H [V 16 5iH] 



(3) 



+ k^ H [V 165 H] 



%^ = ^[V 165 ][H]-k off [V 165 H] 

at K EC m (4) 

-kp[Vi 65 H]-fcHv 165 H] 



Boundary conditions for VEGF reaction at the cell surface 

VEGF and cleaved VEGF are coupled to the cell-surface 
VEGF receptor population through the specification of 
boundary conditions. These boundary conditions equate 
the flux of VEGF perpendicular to the cell surface, J = 
-D-3V/3n, with the rate of reaction at the cell surface. 
For VEGF 165 and VEGF 189 , these reactions include bind- 
ing to VEGFR2, VEGFR1, and NRP1, while for VEGF 114 
and VEGF 12 i, binding occurs to VEGFR2, VEGFR1, and 
the VEGFR1-NRP1 complex (rate parameters for the 
individual isoforms are given in Table 2). 



D 



165" 



a[v 



165 J 



dn 



= kI n R2 [Vi65][R2]-k^ ff R2 [V 165 R2] 

.k^lV^llRll-k^lV^Rl] (5) 
+ C N1 [V 165 ][Nl]-k^ ff N1 [V 165 Nl] 



„ 3[Vll4] , V , R2 
Un 4 = K. 



dn -4 n R2 [V 114 ][R2]-k^ R2 [V 114 R2] 

+ k: n R1 [V 114 ][Rl]-k^ R1 [Vn4Rl] (6) 
+k: n N1 [Vn 4 ][RlNl] - k^ N1 [V 114 RlNl] 

Equations describing the VEGF receptors are given as 
a system of ordinary differential equations over the 
sprout surface. In addition to describing VEGF binding 
to the receptors, they describe the continual balance 
between secretion of unligated VEGF receptors (at a 
rate s R ) and internalization (k int ), VEGF 165 bridging of 
VEGFR2 and NRP1, VEGFR1/NRP1 coupling (k c ) and 
uncoupling (k uc ), and VEGF 114 binding to the VEGFR1- 
NRP1 complex. Note that VEGF 121 and VEGF 189 terms 



Table 2 Kinetic Parameters 



Isoforms Involved 


Reaction 


Forward 

Kin/ K:/ kint 


Reverse 

koff/ kuc 


Keq 


Ref. 




Cell-surface Reactions 










121, 165 


Binding VEGF to VEGFR2 


MO 7 M~Y 1 


0.001 s" 1 


100 pM 


[42] 


All 


Binding VEGF to VEGFR1 


3-1 0 7 M~Y 1 


0.001 s" 1 


33.3 pM 


[42] 


165, 189 


Binding VEGF to NRP1 


3.2-1 o 6 (vrY 1 


0.001 s" 1 


312 pM 


[42] 


165 


Binding VEGF 165 -VEGFR2 to NRP1 


3.1-10 6 
[10 15 um 2 /(mol-s)] 


0.001 s" 1 


323 

mol/10 15 Mm 2 


[42] 


165 


Coupling VEGF 165 -NRP1 to VEGFR2 


MO 7 
[10 15 um 2 /(mol-s)] 


0.001 s" 1 


100 

mol/10 15 |jm 2 


[42] 


121, Cleaved 


Coupling VEGFR1 to NRP1 or 
Coupling VEGF 114 -VEGFR1 to NRP1 

Internalization, k jnt 

Insertion, s R2 , s R1 , or s N1 

ECM Reactions 


MO 7 
[10 15 Mm 2 /(mol-s)] 

2.8-1 0" 4 s" 1 

e.g., k int .-[R2] total 


0.01 s" 1 


1000 
mol/10 15 Mm 2 


[42] 

[42] 
calculated 


165 


VEGF 165 to HSPG 


6.06-1 0 4 M~Y 1 


0.01 s" 1 


165 nM 


[63] 


189 


VEGF! 89 to HSPG 


1.18-10 6 (vrY 1 


0.01 s" 1 


8.5 nM 


Figure S1.1 


All 


Degradation rate, k deg (only Figs. 4C,D; 8A) 


0.001 s" 1 






assumed 


All 


Protease activity, k P (only Figs. 6; 7A; 8A) 


2.8-1 0" 4 s" 1 






assumed 
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are omitted for clarity. We assume that receptors do not 
diffuse along the cell surface. 

d[Rl] c ,Ri mil 
-Sri -k int [Rl] 



dt 



165,R1 



-^[V 165 ][Rl] + k^ Ki [V 16 5Rl] 

&BM (7) 



114,R1 



-[Vii 4 ][Rl]+k^ R1 [V 114 Rl] 
k^^RlKNll+^'^tRlNl] 



Kbm 

R1,N1 r 



^[V 165 R2N1] 
dt 



-k 1 



1 [V 165 R2N1] 



+k 



'[v 165 R2][Ni]-k;; 



165R2.N1 



[V 165 R2N1 



(14) 



+k 165Nl,R2 [Vi65N1][R2] _ kl«Nl,R2 [Vi65R2N1] 



d[V 114 Rl] 



l [V 114 Rl] 



k 114, Rl 

+ ^^[Vn4][Rl]-k^ R1 [V 114 Rl] 



(15) 



Kbm 



dim „ ,R2 rpol 
— =s K -yR2] 



i 165,R2 
K on 

Kbm 

i 114,R2 
K on 

Kbm 



[Vxssl^l + k^^IVifisRl] 



[V 114 ][R2] + k^ R 2[ Vll4 R2] 



(8) 



_ k 165Nl,R2 [Vi65N1][R2] + k J«Nl,R2 [Vi65R2N1] 



i 165,N1 

[V 1S 5][N1] + C 5 ' N1 [V 165 N1] 

- kJ S5R2 ' N1 [V 165 R2][Nl] + kif R2 ' N1 [V 165 R2N1] 

- kf 4R1 ' N1 [V 114 R1][N1] + C R1 ' N1 [ViwRlNl] 
-k R1 ' N1 [Rl][Nl] +k R c 1 ' N1 [RlNl] 



d[R1N1] =-C-[RiNi] 



dt 



fk K1N1 [Rl][Nl]-C N [R1N1] 

^.114, R1N1 

-[V 114 ][R1N1] + k^ ff 4 ' R1N1 [V 114 R1N1] 



(10) 



Kbm 



4Vl65R1] =-k?r i [v 165 Rii 



dt 



i 165, Rl 

+ ^^[V 165 ][Rl]-k^ R1 [V 165 Rl] 



^[V 165 R2] ,165R2r v R?1 

- ~ k int LVi 65 R2J 



(11) 



dt 



i 165, R2 

+ ^ [Vi 65 ][R2] - C 5 R2 [V 165 R2] 



(12) 



k 165R2,Nl [ Vl 65 R 2 ] [N 1 ] + C SKiN ' [Vi 65 R2N 1 



165R2,Nli 



d[V 165 Nl] 



dt 



-k^ N1 [Vl6 5 Nl] 



+ ^ [V 165 ][N1] -k^' N1 [V 165 Nl] 

&BM 

_ k 165Nl,R2 [Vi65N1][R2] + kl65Nl,R2 [Vi65R2N1] 



(13) 



d[Vn 4 R2] i.ii4R2 rv R91 

= - k int l V H4R^J 



dt 



i 114, R2 

+ [V 114 ][R2] - k^ 4 ' R2 [V U4 R2] 



d[V 114 RlNl] = _ CR1N1[Vii4RiNi] 



(16) 



dt 



+-SL [V 114 ][R1N1] -C' R1N1 [V 114 R1N1] (U > 

&BM 

+k 114Rl,Nl [Vn4R1][N1] _ k 114Rl,Nl [Vn4R1N1] 



External Boundary conditions for VEGF secretion and 
clearance 

The computational domain is assumed to be representa- 
tive of its surrounding tissue in the radial direction. 
Thus, we use a no-flux condition, dV/dr = 0 at r = 
R edge . The VEGF gradient was specified by secretion of 
uncleaved VEGF at z = +L using the Neumann BC, 
-D v -3[V]/3z = -q, and a first-order VEGF clearance at z 
= -L, -D v -3[V]/3z = -kciearlY]. The secretion rate q and 
kciear were pre-calculated to generate the desired VEGF 
distribution in the absence of proteases, i.e. a concentra- 
tion V 0 at z = 0 and a VEGF gradient, g 0 over the 
domain length. For cleaved VEGF, there was no secre- 
tion, but there was clearance. Protease activity and 
degradation were assumed to be uniform. 

The volume-averaged clearance rate of soluble VEGF 
in our model, k c i ear /(2L) ~ 2.49-10" 4 s" 1 , is similar in 
magnitude to clearance times in vivo [51] and represents 
internalization by downstream cells and transvascular 
permeability. 

Steady-state assumption and consequences 

The steady-state approximation assumed in our simula- 
tions is justified because VEGF diffusion and reaction, 
which our model considers, occur much more rapidly 
than the structural changes in vessels and parenchymal 
cells, which may take days to weeks. Transients in the 
VEGF distribution are limited by the residence time of 
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VEGF in tissue, which has been experimentally [51-53] 
and theoretically [42] shown to be ~1 h or less. The 
steady nature of the VEGF distribution has been experi- 
mentally confirmed in the hindbrain, where the VEGF 
distribution is relatively unchanged over an 18 h win- 
dow [13] and still maintains isoform specificity. If transi- 
ent diffusion does contribute to the specificity of VEGF 
gradients, we would expect vascular patterning to lose 
its isoform-dependence over long periods, which is not 
seen either in the retina (1 week [12]) or in tumors (sev- 
eral weeks [6,7,19]). The same may be true for in vitro 
and ex vivo systems (1-3 days, [7,54]). 

The VEGF transport equations can be simplified at 
steady state. We perform the analysis for VEGF 165 , but the 
analysis is applicable to any isoform. In the absence of pro- 
teases, VEGF is in equilibrium with HSPG (k on [V][H] = 
k off [VH]) thus the VEGF diffusion equation reduces to: 



V 2 [V] = 0 



(20) 



This equation states that VEGF transport is not influ- 
enced by HSPG binding. In the presence of proteases, 
VEGF 165 can be converted to VEGF 114 , and summation 
of the VEGF 165 and VEGF 114 equations and neglecting 
minor variations in the intrinsic diffusivity of the two 
isoforms, leads to: 



v 2 ([v 165 ] + [v 114 ]) = o 



(21) 



Thus, the total soluble VEGF is specified by the same 
equation regardless of the presence or distribution of 
VEGF-cleaving proteases. If receptor-mediated internali- 
zation is identical between the isoforms (e.g. in the 
absence of NRP1), total soluble VEGF should not be 
altered by proteases. 

Another important quantity, matrix-bound VEGF, can 
be calculated using the relation, [H] = [H] Total - [VH]. 
Assuming a uniform protease distribution, with proteo- 
lysis rate k P , we obtain 



[VH] 



k on [V][H]. 



iTotal 



koff + k P + kon[V] 



(22) 



Note that for typical parameters used in this study, k off 
>> k P >> k on [V], we can approximate VEGF/HSPG 
binding by [VH] ~ [V] [H] Total /K d (i.e. that [VH] is 
directly proportional to [V] and total HSPG). The 
steady-state distribution of the soluble fraction of VEGF 
can thus be described by: 



,[H], 



Total 



1 1 r kp / kn n 

— V 2 V = — 1 + , 
[V] 11 D V k off + kp 



_kp 

DEff 



^ (23) 



Thus, VEGF proteolysis is enhanced not only by the 
rate of proteolysis, k P , but also by matrix binding, which 
decreases VEGF's effective diffusivity, D Eff = D/(l + [H] 



Totai/Kd)- As a result, matrix binding potentiates VEGF 
proteolysis. 

Numerical Methods for VEGF Calculations 

The transport of VEGF and the dynamics of HSPG were 
solved using the finite volume method in cylindrical coor- 
dinates (z and r), Figure 2. The control volume spacing in 
the z direction was 8 \im. In the r direction, one voxel was 
used to represent the sprout radius, from r = 0 to r = 
R sprout ; for r > R spr0 ut (outside the sprout surface), spacing 
was 4-8 |im. A finer grid produces consistent results differ- 
ent by less than 5% and does not alter qualitative conclu- 
sions. Heterogeneous reactions at the cell surface were 
approximated as homogenous reactions, using a basement 
membrane layer formulation given previously [40] . 

The steady-state solution was obtained by solving the 
transient solution from the initial conditions until a rela- 
tive convergence of 10" 7 was achieved at each node. First- 
order temporal derivatives were discretized using a first 
order fully-implicit scheme, while the second-order spatial 
derivatives were discretized using a second order central 
difference method. Nonlinear solution of the equations 
was found by iteration using the successive over-relaxation 
(SOR) update formulation and a Red-Black node ordering 
[55]. Simulations were run on a personal computer using 
custom code written for Matlab 7.6.0. 

Model Implementation and Initial Conditions 

Because the cellular internalization of VEGF creates a 
region of VEGF depletion near the sprout, to be consis- 
tent we ran each simulation with an initial VEGF gradi- 
ent in the absence of proteases or degradation. We first 
imposed Dirchelet BCs satisfying the overall domain 
gradient of VEGF, g 0 , and mean VEGF concentration V 0 
at z = 0 (see Parameters): at z = +L, [V 165 ] = V 0 -(l 
+g 0 -L/L tip ), at z = -L, [V 165 ] = V 0 -(l-g 0 -L/L tip )), setting 
any degradation and proteolysis rates to be zero. This 
simulation can be done in the absence of HSPGs to 
reach convergence more quickly since HSPGs have no 
influence on steady-state soluble VEGF under these con- 
ditions (see previous section). After equilibration of dif- 
fusion and the VEGF receptors, Dirichlet BCs were 
converted into Neumann BCs specified by q and k c iear> 
defined using the following formulas: 



R edge R edge 

q- j rdr= J D 165 



9[V 



165 J 



dz 



rdr 



z=+L 



Redge 



kele 



/ 



[Vi 65 ]rdr = 



Redge 



a[v 



165 J 



dz 



(18) 



rdr (19) 



z=-L 
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Due to averaging over the respective boundaries, these 
equations only approximately specify the conditions, g 0 
and V 0 ; however this is sufficient for our study. After 
this step, we incorporated HSPGs and VEGF-HSPG 
complexes at equilibrium, any degradation terms and 
proteolytic reactions to arrive at the predicted steady- 
state spatial distribution. 

Geometrical and Transport Parameters 

The tissue (cylindrical computational domain) is 160 \im 
long (L = 80 |im) and 100 |im in diameter (R e d ge = 50 
(im) (Table 3, Figure 2A). The sprout is a cylinder from 
z = -80 (im to z = 0 with a radius of R spr out = 2 [im. 
These dimensions were chosen to reflect the average 
sprout size and sprout-to-sprout distance based on 
micrographs of sprouting vasculature [12,13,56]. 

The matrix composition affects both the diffusivity 
and the available volume for VEGF transport. We esti- 
mate the available volume fraction of the ECM as 0.85 
consistent with experimental data on transport studies 
of dextran in cellularized agarose gels [57]. Diffusivity 
was calculated by estimation of the aqueous diffusivity, 
D aq , at 37°C using data from Berk et al. and the Stokes- 
Einstein relation [58,59]. The hindrance of the matrix 
was then approximated using Ogston's relation, assum- 
ing a matrix composition of 14% v/v collagen fibrils (r h 
= 20 nm), 0.078% v/v GAG chains (r h = 0.55 nm) as in 
[60]. Finally, we corrected for the increased viscosity 
due to proteinacious interstitial solution (20.6 g protein/ 
L) [61], resulting in D 165 = 68.6 |im 2 /s. Despite a slight 
error, we standardized the results of the different VEGF 
isoforms by assuming D 165 for also VEGF 114 , VEGF 12 i, 
and VEGF 189 . Receptor movement is neglected as diffu- 
sivity along the cell surface, D ~ 0.01 (im 2 /s [62] is 
much slower than reactions at the cell surface; the Dam- 
kohler number is Da = k 0 ffL tip 2 /D = 160. 

Kinetic Parameters for Reactions 

To estimate the binding of VEGF 165 to HSPGs, we used 
binding rates of bFGF to HSPGs (k on = 4.2-10 5 M'Y\ 
k 0 ff = 0.01 s~\ K d = 23.8 nM) [63] and adjusted k on to 
match K d estimates of VEGF 165 binding to heparin, K d 
~ 165 nM [38,64]. VEGF 189 binds to heparin and 



HSPGs more strongly than does VEGF 165 , however no 
estimate for its K d had been made. Based on extrapola- 
tion of VEGF 165 heparin elution data (described in more 
detail in Additional file 1, section SI; Figure SI. IB), we 
estimated a lower limit for K d of 8.5 nM. VEGF 12 i and 
VEGF 114 do not bind to HSPG [23]. The overall rate of 
VEGF cleavage by proteases in tissues is not known; but 
we assumed a rate that produces significant effects for 
our model parameters, k P = 2.8-10" 4 s" 1 (time-scale of 
1 hr), which is similar to VEGF clearance rates in vivo. 

We previously characterized kinetic parameters for the 
binding of VEGF to VEGFR1, VEGFR2, and NRP1, and 
for the coupling of the receptors (see Table 2) [42]. 
While VEGF 165 has several fold greater affinity for 
VEGFR2 (in the presence of HSPG) than VEGF 121 
[23,38,65] and VEGF 121 may bind NRP1 with low affi- 
nity [45], we assume that the differences in the isoform 
binding to receptors are due to differences in NRP1 
binding: no binding for VEGF 12 i; affinity for VEGF 189 
higher than that of VEGF 165 in proportion to the iso- 
forms' HSPG affinity. The behavior of cleaved VEGF 110 
(and presumably VEGF 114 ) has been shown to be identi- 
cal to that of VEGF 12 i [23] indicating similar receptor 
kinetics to VEGF 121 . Finally, we assume receptors are 
internalized at a rate, k int , independent of VEGF binding 
[66] and that each VEGFR is maintained at a constant 
level, i.e. s R2 = k int -[R2] Total . 

As described above, we studied a hypothetical conti- 
nuum of VEGF isoforms with properties intermediate 
between the three isoforms: we varied VEGF receptor 
binding parameters along this continuum proportionally 
to the isoforms' affinity for HSPGs. Two opposing effects 
needed to be accounted for: the ability of VEGF to couple 
VEGFR2 and NRP1 and the ability of VEGF to bind to the 
VEGFR1-NRP1 complex. Thus the VEGF-NRP1 binding 
rate was k on v ' N1 = k on 165 ' N1 .(K d 165 ' H /K d v ' H ). The coupling 
rate, k c , between VEGF-VEGFR2 and NRP1 was also 
increased by the same factor to preserve equilibrium. In 
contrast, VEGF binding to the VEGFR1-NRP1 complex 
and NRP1 coupling to the VEGF-VEGFR1 complex was 
higher for shorter isoforms, with k'forward = kforward*(Kd V ' H 
- K d 165 ' H )/K d 165 ' H and k reverS e unchanged, for isoform 
HSPG binding affinities greater than K d 165,H = 165 nM 



Table 3 Physical and Transport Parameters 



Parameter 


Value 


Ref. 


Radius of domain, R 


50 Mm 


Assumed 


Half-Length of domain, L 


80 Mm 


Assumed 


Length of tip cell, L cell 


40 Mm 


[12] 


Radius of sprout, R spr out 


2 Mm 


[12] 




Tip cell surface area: 515 Mm 2 


calculated 


Diffusivities, D 1H D 12 i, D 165 , D 189 


68.6 Mm 2 /s (see Methods) 


estimated 


Local Available Volume Fraction of the extracellular matrix, K ECM 


0.85 


[57] 
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(any isoform with greater affinity to HSPGs than VEGF 165 
has no binding to VEGFR1-NRP1). 

A degradation rate constant of k deg = 10" 3 s" 1 was used 
to estimate the impact of significant degradation. As we 
will see, total degradation is independent of the isoform 
if the HSPG-bound isoforms are protected against degra- 
dation. Conversely, degradation is isoform-specific when 
the heavier isoforms can be degraded when bound to 
HSPGs as well as in solution. The heavier isoforms will 
be degraded more due to their longer residence time. 

Concentrations of Receptors and HSPG 

Concentrations of VEGF, abluminal VEGF receptors, 
and ECM HSPG were estimated previously [8,42,60,67] 
(see Table 4). VEGF levels were taken to be V 0 = 1 pM 
as previously determined from the literature [67]. We 
assume that each endothelial cell has a fixed population 
of VEGFR maintained by a balance of receptor insertion 
into the plasma membrane and internalization. We 
assumed a base condition of [Rl] Total = 10 4 an d [R2] Total 
= 10 4 per vascular surface corresponding to the area of 
a tip cell (515 |im 2 ). To account for the excess of NRP1, 
we use [Nl] Total = 3-10 4 molecules/515 \im 2 . The bind- 
ing site density of HSPG in the ECM was 750 nM [41]. 

Definition of Gradient Metrics for VEGF isoforms and 
VEGF receptors 

Endothelial cells may sense absolute or relative gradients 
of VEGF [68,69]. For initiation of vascular sprouting, the 
gradient across one cell is likely to be sensed. For sprout 



extension and maintenance of the tip cell/stalk cell dif- 
ferentiation, the signaling difference between two cells 
(e.g. a tip and stalk cell [70]) may be more relevant. 

In the present study we consider gradients only 
between the front and back of the surface of the tip cell, 
with mean concentrations ([V] T i P ) defined over the tip 
cell's entire surface. The absolute gradient, denoted by 
AG V or A T i P [V] is the difference in VEGF concentration 
between the tip cell front and back. Note that the abso- 
lute gradient is additive over all isoforms, e.g. AG V = 
AG165 + AG121 for a two isoform distribution. The 
steepness of the VEGF distribution (relative gradient) is 
the absolute gradient normalized to either the mean tip 
total soluble VEGF, i.e. fractional gradient or FGi 65 , or 
to the mean tip concentration of the isoform under con- 
sideration, e.g. IFG165: 



IFG16 



AGi 



65 



[Vi65] T ip 
1 

[Vi65] T ip 



(24) 



A T ip[Vi 65 ] « A Tip ln [Vi 65 ] 



FG 



AG 16 



165 



[Vi 65 ] Tip + [Vi2l] Tip 
1 



(25) 



[V 165 ] Tip + [V 121 ] 



-A Tip [V 



165 J 



Tip 



Note that the above metrics differ only in the presence 
of multiple VEGF isoforms. 



Table 4 Model Conditions 



Parameter 


Symbol 


Value 


Ref. 


Initial mean VEGF level 


V 0 


1 pM 


[67] 


VEGF gradient over domain 


g 0 


5%/40 Mm 


Assumed 


[HSPG] in ECM 


[H] Total 


0.75 -10~ 6 mol/L 


[41] 


System with receptors 








Total VEGFR2 levels per area of tip cell on vasculature 


[R2] Tota i 


10,000/(area of tip cell); area of tip cell = 515 urn 2 


[42] 


Total VEGFR1 levels per area of tip cell on vasculature 




10,000/(area of tip cell) 


Assumed 


Total NRP1 levels per area of tip cell on vasculature 


[NUrotal 


30,000/(area of tip cell) 


Assumed 


VEGF (typically VEGF 165 ) secretion rate 


q 


5.27-1 0" 5 molec/um 2 -s 


Calculated 


VEGF clearance rate 


Kdear 


0.0399 um/s 


Calculated 


System without receptors 








VEGF (typically VEGF 165 ) secretion rate 


q 


4.39-1 0" 5 molec/um 2 -s 


Calculated 


VEGF clearance rate 


Kdear 


0.0948 um/s 


Calculated 


HSPG gradient (Figure 4B) 








HSPG gradient 


9ho 


30%/40 Mm 


Assumed 


Soluble VEGF degradation (Figure 4C) 








Degradation rate of soluble VEGF 




0.001 s" 1 


Assumed 


Degradation rate of matrix-bound VEGF 


K deg 


Os" 1 


Assumed 


Matrix-sequestered VEGF degradation scheme (Figure 4D) 








Degradation rate of soluble VEGF 


kdeg 


0.001 s" 1 


Assumed 


Degradation rate of matrix-bound VEGF 


Meg 


0.001 s" 1 


Assumed 
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We use similar definitions to find gradients of cell- 
surface VEGF-bound receptor complexes. In this case, 
the gradient is limited to the axial (z) direction due to 
rotational symmetry in our model. We normalize the 
absolute gradient of receptor complexes to the total 
amount of a particular receptor, e.g.: 



AGi 



1 



65-R2 



[R2J 



Total, Tip 



-A T ip[V 165 R2] 



(26) 



where [R2] Total> Tip is the sum of all forms of VEGFR2 on 
the tip cell surface. Using the fractional occupancy (FO) of 
a VEGF isoform to a VEGF receptor, the fractional gradi- 
ent of VEGF 165 bound to VEGFR2 in a system consisting 
of isoforms VEGF m and VEGF 165 can be expressed as: 



FGi 



65-R2 = 



[V 165 R2] Tip + [V 121 R2] Tip 

AG165 - R2 



A Tip [V 165 R2] 



(27) 



roi 6 5_j?2 + FO12I- 



R2 



We introduce the concept of sensitivity, S, to define 
the multiplicative relation between a metric in VEGFR2 
binding and a similar metric in the VEGF distribution, 
e.g. S IFG _ 165 = IFG 165 _ R2 /IFG 165 . Note that S IFG also 
describes the factor increase in the fractional occupancy 
(FO) given a factor increase in [VEGF], strictly valid 
only in the differential limit: 



SlFG-l 



/ 1 d[V 165 i?2] \ / 1 
V[V 165 i?2] dz )'\[Vi6i 
( d[V 165 R2} \ ( d[V 165 ] y l d_ 
\ \V 165 R2] ) ' \ [V 165 ] ) d 



dz 
ln(FO) 
ln[Vi 65 ] 



(28) 



For a single receptor detecting a single ligand, S IFG is 
always less than 1 (S IFG = 1 - FO) (f° r more detail, see 
Additional file 1, section S2.2). 

Additional Assumptions 

We assume that degradation and proteolysis is uni- 
formly distributed. Our model does not also fully reflect 
the contribution of vasculature to VEGF transport. 
While vasculature-mediated internalization is thought to 
represent the primary clearance for VEGF in tissue [43], 
our model only implicitly represents such effects 
(through a clearance at z = -L), neglecting isoform-spe- 
cific variation. We also assume constant internalization 
and insertion rates for the VEGF receptors, resulting in 
constant total receptor population on the sprout. This is 
justified because we are primarily interested in the sen- 
sory role of the sprout and not in its internalization 
capacity, however ligand-induced shifts in NRP1 and 
VEGFR2 may also alter sensory function. Mechanistic 
models of these processes have not been developed and 
thus we did not incorporate them in the current study. 



Results 

The computational model predicts that the steady-state 
distribution of soluble VEGF is isoform-independent when 
considering only diffusion and matrix binding 

When VEGF is secreted, extracellular HSPGs are pre- 
dicted not to influence the distribution of the soluble 
fraction of VEGF (Figure 3A vs. 3B); nor do they influ- 
ence receptor signaling (Figure 3, bar graphs), assuming 
VEGF-HSPG complexes in the ECM cannot directly 
ligate VEGF receptors. At steady state, each molecule of 
VEGF that binds to the matrix is matched by a molecule 
that unbinds. Note that using physiological parameters, 
VEGF concentration is significantly below the K d for 
matrix binding ([V] ~ 1 pM, K d = 165 nM) resulting in 
HSPG binding sites remaining unsaturated; in this 
regime, matrix-bound VEGF concentration is linearly 
proportional to soluble (unbound) VEGF concentration 
at all spatial positions. 

At the tip cell, the absolute gradient of unbound 
VEGF is -0.059 pM/40 urn. The direction of the gradi- 
ent of receptor signaling (characterized in this study by 
the receptor occupancy) is the same, with magnitude 
-30 VEGF-VEGFR2/40 urn. The relative steepness (iso- 
form fractional gradient, IFG) of both the VEGF and 
VEGF-VEGFR2 distributions are similar, IFG V egf = 
6.39%/40 urn and IFG V egf-vegfr2 = 6.1296/40 urn. Thus, 
the sensitivity of VEGFR2 binding to VEGF gradients 
(Sifg) is high (Figure 3B). Note that VEGF receptors 
also locally deplete VEGF and thereby influence the dif- 
fusion of VEGF (Figure 3, arrowed dotted lines); in the 
case shown, receptors capture -39% of the secreted 
VEGF. 

Differences in VEGF gradients can arise from differences 
in VEGF degradation 

The lack of influence of HSPGs on VEGF gradients at 
steady state suggests that different isoforms secreted at the 
same rate (as occurs in transgenic systems designed to 
express only one isoform), and differing only in HSPG 
binding, will have similar soluble VEGF distributions 
(Figure 4Aii). Matrix-bound VEGF increases with HSPG 
affinity of the isoform (Figure 4Aiii), as does the total 
amount of extracellular VEGF (for VEGF 12 i: 646 mole- 
cules; VEGF 165 : 3.58-10 3 molecules; VEGF 189 : 5.76-10 4 
molecules). However, because bound VEGF is propor- 
tional to soluble VEGF, the relative gradients of matrix- 
bound VEGF are the same for each isoform. Over the 
simulated length scales, total VEGF concentration gradi- 
ents of different isoforms would not intersect (Figure 
4Aiv), regardless of HSPG concentration {not shown). 
Thus, differential HSPG binding of isoforms is not suffi- 
cient to account for observed VEGF patterning in vivo 
(Figure 1A). 
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Figure 4 Mechanisms of forming differential isoform gradients. We consider the distributions of VEGF 12 i, VEGF 165 , and VEGF 189 under four 
scenarios that may be responsible for isoform patterning in vivo (i, soluble fraction; ii, bound fraction; iii, soluble + bound VEGF). A, reversible 
HSPG binding considered previously in Figure 3 (HSPG-binding-only model). B, patterning of the underlying HSPG ([H] Tota i = 750 nM at z = 0, 
30%/40 urn). C, soluble VEGF degradation; this is isoform-independent degradation (all isoforms are degraded at the same rate), but HSPG- 
bound VEGF is protected from degradation; k deg = 10~ 3 s" 1 . D, matrix-sequestered VEGF degradation; all isoforms are degraded, and HSPG 
binding confers no protection; this has the effect of increasing degradation of HSPG-binding isoforms due to their longer residence time, even 
though all isoforms have the same degradation rate constant k deg = 10" 3 s" 1 . 



We simulated whether a patterned ECM (Figure 4B) 
could cause the observed gradients, but this scenario also 
did not lead to isoform-specific differences in soluble 
VEGF; however, it did lead to larger isoforms having 
steeper distributions in bound and total VEGF (Figure 
4Biv). Still, VEGF levels between different isoforms are 
not predicted to intersect within the spatial frame of 
reference. 

We first considered whether the steady-state assump- 
tion was accurate, since transient diffusion of isoforms 
having different HSPG affinities can momentarily lead 
to differential isoform patterning. However, experimen- 
tal evidence points to VEGF 12 i and VEGF 165 having 
rapid kinetics in vivo (x ~ 1 h) [51,71], which likely indi- 
cates that VEGF patterning operates close to steady 
state. We will consider the ability of VEGF-cleaving pro- 
teases to create these spatial distributions; however, it is 
of note that VEGF fragments have not yet been 
observed in developing tissues [17]. 



Next, we considered VEGF degradation (Figure 4C,D). 
In vivo, VEGF diffuses in an environment filled with 
cells, which may selectively uptake VEGF isoforms, e.g. 
through cell surface receptors, possibly in conjunction 
with NRP1 or cell-surface HS [72]: these effects are 
incorporated into the VEGF degradation term of our 
model. Degradation may also result from VEGF inactiva- 
tion by isoform-selective VEGF inhibitors, e.g. connective 
tissue growth factor, sVEGFRl, thrombospondin-1, or by 
proteases that can cleave VEGF to fragments that are not 
recognized by commonly employed antibodies. 

As described in the methods section, we considered 
two specific cases for degradation: first, the "soluble 
VEGF degradation moder in which HSPG-bound VEGF 
is protected from degradation (Figure 4C); and the 
"matrix- sequestered VEGF degradation" model, in which 
HSPG-bound VEGF can be degraded (Figure 4D). The 
former is similar to heparin-bound bFGF being pro- 
tected against plasmin cleavage [73]. 
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Computational simulations suggest that when HSPGs 
protect VEGF (Figure 4C), the soluble VEGF distribu- 
tion is independent of HSPG binding, as if the HSPGs 
were not even present. The relative gradients of the 
VEGF isoforms are equal (8.4%/tip length) and thus can- 
not recapitulate isoform-dependent localization. In addi- 
tion, total VEGF in the system is approximately 
proportional to the HSPG binding affinity. 

In contrast, when HSPG-bound VEGF is vulnerable to 
degradation (Figure 4D), the different grading of VEGF 
isoforms seen in vivo (e.g. Figure 1A) is recapitulated 
[7,21,22,39]. Specifically, the strongly heparin-binding 
VEGF 189 produces a very steep distribution (Figure 4Div). 
As a consequence of this, close to the source of secretion, 
increased ECM-binding affinity (i.e. longer isoform) corre- 
sponds to increased levels of bound VEGF; but at suffi- 
ciently far distances, this trend is reversed and the highly 
diffusible isoform is more abundant (Figure 4Diii). Despite 
different spatial localization, the total amount of VEGF in 
the system is similar for each isoform (for VEGF 12 i: 229 
molecules; VEGF 165 : 345 molecules; VEGF 189 : 379 mole- 
cules), and may be consistent with data from [7,22]. 

Proteolytic cleavage of bound VEGF increases soluble 
VEGF only in a model of isoform-specific degradation 

VEGF cleavage (generating sizeable fragments) is also an 
important modulator of vascular patterning and experi- 
mentally, increased cleavage seems to mimic the effects of 
reduced HSPG binding affinity: increase in soluble VEGF, 
increase in range of the VEGF signal with a simultaneous 
decrease in the magnitude of the VEGF gradient, decrease 
in matrix-bound VEGF, and possibly, preservation of the 
total tissue VEGF levels [7,12,13,17,19,21]. Here we con- 
sider the role of VEGF-cleaving proteases on the VEGF 
distribution of the different isoforms for two different sce- 
narios: the HSPG-binding-only model, which provides a 
reference case (Figure 5) and the model including matrix- 
sequestered VEGF degradation (Figure 6), which we saw 
above can result in VEGF isoform gradients similar to 
those observed experimentally. We assume VEGF cleavage 
can act on both soluble and matrix-bound VEGF. 

Surprisingly, inclusion of proteolytic cleavage did not 
affect total soluble VEGF levels in the HSPG-binding- 
only model, despite reducing levels of uncleaved (full- 
length) VEGF and producing VEGF 114 , a non-HSPG 
binding diffusible isoform (Figure 5Ci,Di,Ei,). Mathema- 
tical analysis indicates that this effect is similar to the 
HSPG binding independence noted earlier. Proteases are 
predicted to reduce matrix-bound VEGF levels (Figure 
5Dii,Eii), due to a proportional decrease in uncleaved 
VEGF, which is accompanied by a significant steepening 
of these distributions (fractional gradients are noted 
next to the lines on Figure 5). Comparing the different 
VEGF isoforms, proteases cleave VEGF 189 much more 



thoroughly than VEGF 165 , while total VEGF in the 
VEGF 12 i system (Figure 5Cii) is not affected by proteo- 
lysis because its cleavage produces an equivalent iso- 
form. This potentiation of VEGF proteolytic cleavage by 
HSPG binding is due to increased residence time of the 
matrix-binding isoforms. 

In contrast, with the inclusion of matrix-sequestered 
VEGF degradation (Figure 6), proteolytic release of 
VEGF resulted in increased soluble VEGF levels, size of 
this increase was proportional to the matrix affinity of 
the uncleaved isoform (Figure 6Ci,Di,Ei). Furthermore, 
both the soluble VEGF and total VEGF distributions 
have reduced gradients (Figure 5Hi,ii) indicating that 
VEGF gradients can be disrupted by proteases in this 
model. Matrix-bound VEGF decreased due to the action 
of protease; however, its decrease was not as significant 
as in the HSPG-binding-only model, due to the exis- 
tence of VEGF degradation processes. Overall, the total 
VEGF amount in each tissue is well preserved by the 
action of proteases (in the presence of proteases, 
VEGF 12 i system: 229 molecules; VEGF 165 system: 380 
molecules; VEGF 189 system: 428 molecules; compare to 
Figure 4D levels; refer to Figure S3.1). 

We note that in each case, cleaved VEGF has a nearly 
flat concentration profile, a consequence of it not being 
secreted at a single location. This contrasts with the 
uncleaved soluble and bound fractions, which are pro- 
portional to each other and are sharpened by proteases. 
Another surprising finding is that in the presence of iso- 
form-specific degradation, the concentration of cleaved 
VEGF is similar in each isoform system (Figure 6Ci,Di, 
Ei). As described in more detail in Additional file 1, sec- 
tion S2, this is a result of total VEGF being preserved, 
indicating that the residence time of each isoform, and 
hence the time available for cleavage to occur, is similar. 

Correlating VEGF patterning to vascular patterning 
suggests total soluble VEGF may be an optimal VEGF 
signal 

Using the VEGF patterning in vivo as recapitulated in the 
model, we now ask whether there is a viable VEGF-based 
guidance cue for cell sensing that is consistent with the 
twin observed properties: monotonic increase in vascular 
density with isoform length; and the opposite effect of 
VEGF cleavage (Figure 7A). We present results for VEGF 
near the tip cell (z = -40 to 0 (im), with an extended 
range of VEGF isoforms and MMP activity levels to com- 
pare, side-by-side, the isoform dependence and the 
MMP-based proteolytic dependence (Figure 7B,C). 

In the HSPG-binding-only model (Figure 7B), total 
soluble VEGF is constant, and thus is unlikely to carry 
relevant information regarding vessel patterning. Instead, 
such information could come from changes in: soluble 
uncleaved VEGF; cleaved VEGF; or matrix-bound 
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Figure 5 Effect of VEGF-cleaving MMPs on the VEGF distribution. Proteases (MMPs) were added to the HSPG-binding-only system at k P = 
2.8-1 0" 4 s" 1 for systems expressing VEGF 12 i (C), VEGF 165 (D), or VEGF 189 (E). We look at the components of the soluble VEGF distribution (i), and 
bound and total VEGF (ii), both before the addition of proteases (solid lines), and after (dashed lines). The total VEGF distributions of different 
isoforms, in the presence of proteases, are displayed against each other (B). Simulations were performed in the absence of receptors (parameters 
given in Table 4). Isoform fractional gradient values at the tip cell are noted in grey with units %/40 urn. 



VEGF. In order for a signal to recapitulate the isoform 
and MMP dependence shown in Figure IB, it should 
behave in the same fashion to decreasing HSPG affinity 
as it does to increasing MMP activity. In the context of 
Figure 7, that means that the behavior moving to the 
right along the x-axis should be the same on the HSPG 
and MMP graphs. Here however, levels of soluble 
uncleaved VEGF increase monotonically as the HSPG 
affinity of the secreted isoform decreases (due to 
decreased residence time in the ECM), but decrease as 



MMP activity increases. Cleaved VEGF (Figure 6B, red 
lines) is also inconsistent. Matrix-bound VEGF, however, 
has a biphasic behavior; at the tip cell, isoforms with 
intermediate and lower affinity to HSPGs satisfy the 
requirement that MMP effects oppose HSPG binding 
effects, but above an affinity limit this is violated (Figure 
6B, yellow shaded region). This biphasic effect arises 
because greater matrix binding affinity results in greater 
accumulation close to the source of secretion, but also 
more rapid deterioration of the signal in space caused 
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Figure 6 Effect of VEGF-cleaving MMPs on the VEGF distribution. Proteases (MMPs) were added to the matrix-sequestered VEGF 
degradation system at k P = 2.8-1 0" 4 s" 1 for systems expressing VEGF 12 i (C), VEGF 165 (D), or VEGF 189 (E). We look at the components of the soluble 
VEGF distribution (i), and bound and total VEGF (ii), both before the addition of proteases (solid lines), and after (dashed lines). The total VEGF 
distributions of different isoforms, in the presence of proteases (B, black lines), are compared to those in the absence of protease (B, grey lines). 
Simulations were performed in the absence of receptors (parameters given in Table 4). Isoform fractional gradient values at the tip cell are noted 
in grey with units %/40 urn. 



by proteolysis over a longer residence time. This effect 
was also seen in Figure 5B, as total VEGF 189 levels fall 
below those of VEGF 165 at the tip cell. Note that the 
closer the tip cell is to the source, the greater the likeli- 
hood that the monotonic behavior between isoforms is 
maintained, indicating that the distance from the VEGF 
source, all other factors being equal, may be a major 
determinant of isoform-specific vascular patterning. 

In the matrix-sequestered VEGF degradation model 
(Figure 7C), soluble uncleaved VEGF, cleaved VEGF, and 
matrix-bound VEGF follow similar behavior as in the 



above HSPG-binding-only model. However, total soluble 
VEGF is not constant in this model. Interestingly, while 
neither soluble uncleaved nor cleaved VEGF are individu- 
ally viable guidance cues, their sum (total soluble VEGF), 
displays the desired behavior, increasing with decreased 
isoform affinity to the matrix and with increased MMP 
activity, though at sufficiently high affinities, the differ- 
ences at HSPG affinities larger than that of VEGF 165 are 
not large. In addition, this behavior is present over the 
entire range of parameters, indicating that it is robust in 
its ability reproduce vascular patterning. 
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Figure 7 Competing influences of HSPG binding and VEGF cleavage on biological responses. Experimental evidence indicates that HSPG 
binding and VEGF cleavage modulate the same axis of vascular phenotypes (A) [6,7,13,17]. To understand how the VEGF distribution may be 
responsible for vessel patterning, we measured various metrics of the VEGF distribution as they would vary either with variation in the HSPG 
affinity of the secreted isoform (with MMP activity held constant at k P = 2.8-1 0" 4 s"\ dash-dotted lines) or with variation in the background 
proteolytic activity (k P ranged from 10~ 6 s" 1 and 10" 2 s"\ assuming that VEGF 165 is secreted). The HSPG-binding-only model (B) and the matrix- 
sequestered VEGF degradation model (C) were considered. Circle markers indicate respective positions of (from left to right) VEGF 189 , VEGF 165 , 
and VEGF 12 i (K d = «>). Control cases simulated in the absence of MMPs are given in solid lines. Yellow boxes in B and C indicate the range 
where bound VEGF does not display behavior consistent with the experimental data. Simulations were performed in the absence of receptors, 
the presence of which would have a negligible impact on the VEGF distribution, to obtain computational efficiency permitting simulation of a 
broad range of parameters (parameters given in Table 4). 
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In both models (Figure 7B,C), total VEGF and total 
uncleaved VEGF (i.e. soluble uncleaved VEGF + matrix- 
bound VEGF) (not shown), display similar biphasic beha- 
viors as does matrix-bound VEGF and thus do not 
enhance the signal provided by the matrix-bound form. 
The model also predicts the magnitude of VEGF gradi- 
ents present near the tip cell. While fractional gradients 
appear to be significantly poorer candidates for guidance 
cues than are VEGF concentrations, absolute gradients 
seem to be as effective as VEGF concentrations (Addi- 
tional file 1, section S3; Figs. S3. 2, S3. 3). Absolute gradi- 
ents greatly extend the range of isoforms for which 
matrix-bound VEGF in the HSPG-binding-only model 
displays isoform monotonicity (Additional file 1, Figure 
S3.2Biii), however in the matrix-sequestered VEGF 
degradation model, it seems to show little variation with 
respect to changes in MMP activity, indicating weak 
HSPG/MMP antagonism (Additional file 1, Figure S3.3). 

Effects of receptor binding on isoform monotonicity and 
HSPG/MMP antagonism 

Endothelial cells respond not to VEGF gradients directly 
but to VEGF receptor activation. Can VEGF receptors 
accurately sense total soluble VEGF, which we identified 
as a potential cue, or do other extracellular signals 
become relevant? Note that in the absence of NRP1 
expression, VEGFR2 and VEGFR1 binding will be simi- 
lar to the soluble distributions given in Figure 7, because 
uncleaved and cleaved VEGF isoforms bind similarly to 
the VEGF receptors. The presence of NRP1 may com- 
pensate for the inability of longer isoforms to diffuse far 
in tissues. 

We imposed the VEGF distributions calculated in Fig- 
ure 7 onto our model of VEGF receptor interactions 
(refer to Additional file 1, section S4.1). VEGFR2 binding 
was similar for the matrix-sequestered VEGF degrada- 
tion model (Figure 8) and the HSPG-binding-only 
model [not shown). VEGFR2 binding shows greatest 
binding at isoform affinities between VEGF 12 i and 
VEGF 165 , violating isoform monotonicity (Figure 8Ai), a 
result of the isoform-specific potentiation of the under- 
lying uncleaved VEGF by NRP1. At high enough matrix 
affinities this is, however, unable to compensate for the 
low levels of uncleaved VEGF seen for these isoforms. 
We note that uncleaved VEGF 189 (and possibly other 
long isoforms) cannot bind VEGFR2 [46], though its 
cleavage product can; this does not significantly alter 
the shape of the curve (Figure 8Ai), as most VEGF 189 is 
cleaved (Figure 7C) 

We note that when the uncleaved isoform exhibits 
strong HSPG binding, total soluble VEGF is due primar- 
ily to cleaved VEGF, and when the uncleaved isoform 
has minimal ECM binding, then cleaved VEGF is not 
significant (Figure 7C). A similar feature was noted in 



VEGFR1 binding in the presence of NRP1 (Figure 
S4.1A). As a result, we also tested the potential role of 
VEGFR1 signaling in recapitulating total soluble VEGF 
signaling (Figure 8B). Total VEGFR1 binding displays a 
greater range of isoform monotonic behavior than 
VEGFR2 binding, however, it is limited only up to the 
VEGF 165 isoform, beyond which VEGFR1 binding 
remains flat or may actually increase with increased 
matrix-binding affinity {refer to Additional file 1, Figure 
S3.3DU). 

Discussion 

Our study proposes an important reinterpretation of the 
role of HSPGs and VEGF-cleaving proteases in the con- 
trol of VEGF patterning in vivo. Our central proposition 
is that differential isoform degradation (or clearance), 
and not a difference in diffusion arising from HSPG 
binding, controls the spatial localization of VEGF in tis- 
sues (Figure 9A). 

In the developing hindbrain system, VEGF 12 o distri- 
butes in a much more disperse manner (Figure 1A) than 
VEGF 164 [13]. We find that while VEGF 164 may diffuse 
slower in tissues than VEGF 120 > differences in grading 
will only arise when VEGF 164 is also degraded more 
rapidly by the surrounding hindbrain parenchyma, possi- 
bly due to specific uptake by cell surface HSPGs or 
NRP1, preventing it from diffusing far (Figure 9A). We 
propose that this difference in degradation rates between 
isoforms (which we term isoform-specific degradation, 
and is here a result of degradation of matrix-bound 
VEGF) also explains why secretion of non-heparin bind- 
ing isoforms such as VEGF 113 or VEGF 12 o leads to 
greater VEGF levels in the solution phase of tissues than 
secretion of heavier isoforms, given identical secretion 
rates [7,21]. Note that assuming a simpler mechanism of 
reversible HSPG binding in an in vivo system is not able 
to explain this behavior at steady state (Figure 4A). 
Importantly, our mechanism of isoform-specific degrada- 
tion also provides a general explanation for the phenom- 
ena of protease-mediated redistribution of VEGF activity: 
proteases, by either cleaving VEGF [7], cleaving HSPGs 
[15], or by cleaving VEGF inhibitors [32], inhibit the 
degradation of VEGF (Figure 9A), thereby leading to the 
accumulation of VEGF in tissues and an increased range 
of receptor activation on the vasculature. This hypothesis 
is distinct from what is commonly noted in the literature 
as VEGF release, in which VEGF is focally released from 
the matrix by the action of proteases and diffuses to acti- 
vate vasculature. 

The second thrust of our current work is the elucida- 
tion of how VEGF isoforms and VEGF-releasing pro- 
teases might regulate vascular morphology (Figure IB). 
We originally asked the question: which VEGF metrics 
may be responsible for guiding vascular patterning, to 



Vempati et al. BMC Systems Biology 201 1, 5:59 
http://www.biomedcentral.eom/1752-0509/5/59 



Page 1 8 of 24 



recapitulate both isoform monotonicity [6,12,13,25] and 
the protease dependence in vascular patterning [6,7]? 
We identified several possibilities including the levels of 
matrix-bound VEGF, soluble VEGF levels, NRP1 depen- 
dent signaling, and even VEGF gradient directionality. 
Using the results of our model, we can now provide a 
more specific answer to this question. 

For example, the model predicts as expected that 
VEGF 188 -secreting tumors would have the greatest levels 
of peritumoral matrix-bound VEGF [6], however, farther 
away at the vascular front, matrix-bound VEGF levels 
may actually be lower than if VEGF 165 were being 
secreted (Figure 9B). Thus, how a sprout experiences 
the ordering of the isoforms by sensing matrix-bound 
VEGF may not always be monotonic with respect to the 
heparin binding affinity of the VEGF isoform; it will be 
dependent on the distance between the sprout and the 



VEGF source (Figure 9B). (Note that the specific order- 
ing depends on HSPG concentrations, protease levels, 
degradation rates, etc. Given their intrinsic variability in 
biological systems, the conclusion that a sufficiently 
strong matrix-binding isoform will have less matrix- 
bound VEGF is robust.) Instead, the possibility that 
VEGF 189 -secreting systems also have the lowest levels of 
soluble VEGF of the different isoform-expressing sys- 
tems seems to be consistent, even in the presence of 
proteases (Figure 9B). 

An important alternative to the concentration of 
matrix-bound VEGF levels may however be its direction- 
ality, especially as VEGF 189 seems to have the sharpest 
distributions in the absence of proteases (Figure 4Diii). 
While this holds true in the absence of proteases and can 
possibly explain observations in the mouse hindbrain 
[13], this metric does not account for the effect of VEGF- 
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Figure 8 VEGFR2 and VEGFR1 signaling in mediating vascular phenotypes. Using the VEGF distributions calculated by the isoform-specific 
degradation model in Figure 7 we calculated VEGF binding to VEGFR2 (A) and VEGFR1 (B). We assume 10 4 VEGFR2, 10 4 VEGFR1, and 3*1 0 4 NRP1 
per tip cell (assuming no VEGF depletion by the sprout). VEGF isoforms show isoform-specific differences in their ability to bind NRP1 and the 
VEGFR1-NRP1 complex. To calculate the receptor binding parameters for arbitrary isoforms, we imposed two constraints to interpolate/ 
extrapolate around the known VEGF isoforms (refer to Methods, Kinetic Parameters). VEGF/NRP1 affinity was made to recapitulate the lack of 
binding of VEGF 12 i and stronger binding affinity of VEGF 189 [37]. VEGF binding affinity to the VEGFR1-NRP1 complex and NRP1 coupling to the 
VEGF-VEGFR1 complex was set such that VEGF 121 can fully bind to VEGFR1-NRP1 [95], while any isoform with greater HSPG affinity than that of 
VEGF 165 cannot bind. Circle markers indicate (from left to right) VEGF 189 , VEGF 165 , and VEGF 12 i (K d = «>). 
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Figure 9 Summary of VEGF transport results. A, We propose that isoform-specific patterning in tissues results from differences in their rates 
of degradation (i.e. isoform-specific degradation), potentially mediated by HSPG or NRP1 binding. For longer isoforms to have more localized 
VEGF distributions [13], their degradation must be more rapid. We hypothesize that proteases enhance the spatial propagation of VEGF (i.e. VEGF 
redistribution) [16,17,19] by inhibiting this isoform-specific VEGF degradation. B, Role of matrix-binding affinity on VEGF patterning. Our 
computational model shows that HSPG-binding affinity affects soluble and matrix-bound VEGF in different manners. Increased matrix binding 
affinity universally decreases total soluble VEGF, in the case of HSPG-mediated isoform-specific degradation. However, matrix-bound VEGF shows 
increased "localization", a result of two competing behaviors: an increase in concentration due to accumulation close to the source of secretion, 
and a decrease in concentration farther away due to an increased rate of degradation. This causes matrix-bound VEGF to lose its isoform- 
monotonic behavior. The ordering of how isoforms are perceived is dependent upon the distance that sensing occurs from the VEGF source. 

s. J 



cleaving proteases, which further sharpen the distribu- 
tion, instead of behaving in a HSPG-antagonistic fashion 
(Figure S4.3A). VEGF 189 also may have stronger NRP1 
binding than VEGF 165 [37], however our results show 
that soluble uncleaved VEGF 189 levels may be so low (the 
majority will be cleaved or degraded), that total NRP1- 
potentiated VEGFR2 signaling may be weaker than that 
in VEGF 165 -secreting systems. In fact, our model shows a 
surprising result: VEGFR1 binding of total soluble VEGF 
seems to better reconcile both isoform monotonicity and 
the antagonistic relationship between HSPG affinity and 
MMP activity than VEGFR2 binding does (Figure 9). 

Our results may provide insight into the nature of sev- 
eral experimental observations regarding vascular pat- 
terning. A comparison of our results to those of 
Ruhrberg et aTs hindbrain data suggest that while the 
overall VEGF distribution is spatially non-monotonic 
(Figure 1A) [13], the underlying soluble fraction of 
VEGF is monotonic and is the basis for endothelial 
behavior. We note, however, that several important bio- 
logical effects need to first be taken into account, for 
example the roles of filopodia [12,13,74] and the direct 
receptor signaling of matrix-bound VEGF [35]. 

The loss of monotonicity for matrix-bound VEGF 
(Figure 7B, Figure 9B) is a result of the isoform-depen- 
dent decrease of the soluble uncleaved VEGF fraction in 
space. We now discuss how the tip cell filopodia can 
contribute to the process of isoform sensing. By project- 
ing out in front of the tip cell, filopodia may be able to 
detect a region of space where matrix-bound VEGF may 
in fact operate in an isoform-monotonic manner (Figure 
9B), effectively increasing the spatial range of the 
matrix-bound VEGF fraction's isoform monotonicity. 



We conceptualize this to sprouting angiogenesis. In the 
initial stages of sprouting where the sprout is far away 
from the VEGF source, only soluble VEGF will exhibit 
isoform monotonicity. However, as a sprout continues 
to invade closer into the VEGF secreting tissue, matrix- 
bound VEGF signaling and VEGF gradients will also 
exhibit the correct isoform monotonicity (refer to Addi- 
tional file 1, section S3.2, Figure S4.4). These effects 
would be relevant in vivo as VEGF gradients never seem 
to be more than 50 urn in front of the vascular front 
[12,13], a range that can easily be sensed by filopodia 
[12]. Thus, while total soluble VEGF is theoretically the 
most isoform-monotonic signal, there may be no practi- 
cal differences between it and matrix-bound VEGF. In 
fact, the information provided in matrix-bound VEGF 
may be even more relevant than that of soluble VEGF 
as it can result in differential VEGFR2 signaling that 
favors activation of p38/MAPK compared to soluble 
VEGF [35] and may mediate increased interaction with 
cell-surface NRP1 [34], both of which may play a direct 
role in branching and migration behaviors leading to 
increased vascular density. Unlike the behavior of 
matrix-bound VEGF, an isoform monotonic signal in 
NRPl-mediated VEGFR2 signaling does not necessarily 
emerge as a sprout invades into the VEGF secreting tis- 
sue (a result of reaction limitations in VEGF/VEGFR2/ 
NRP1 coupling) (Figure S4.4C), indicating that NRPl- 
mediated VEGFR2 signaling may be less capable of giv- 
ing rise to isoform-specific increases in vascular branch- 
ing complexity. However, this conclusion is based on 
assuming a model where receptors only bind soluble 
VEGF; how the signaling of matrix-bound VEGF 
changes these conclusions is not known. 
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Our results indicated that VEGFR1 signaling of soluble 
VEGF may be able to also mediate isoform-specific dif- 
ferences in vascular patterning. However, this conclusion 
is at odds with the VEGFR2-dependent nature of angio- 
genesis in the retina [12], in in vitro patterning of porcine 
aortic endothelial cells [7], and in MMP9 induced carci- 
nogenesis [16]. In addition, VEGFR1 signaling supports 
migratory behavior through p38/MAPK [26], however 
our results suggest an inhibitory role of VEGFR1 signal- 
ing on vascular density (Figure 8Bi). While our model 
cannot elucidate the importance of VEGFR1 signaling, it 
is interesting to note that when VEGFR2 signaling is 
monotonically increasing with the isoform length, 
VEGFR1 signaling is decreasing, suggesting that a bal- 
ance between the two receptors may be important. 

The unique role of matrix-bound VEGF in mediating 
the branching phenotype through filopodia may offer an 
interesting solution to the paradoxes of VEGF-cleaving 
MMPs in tumor growth. For example, both VEGF 188 
tumors and VEGF 164A108 _ 118 tumors are similar in that 
they have higher intratumoral vascular density and 
lower-diameter vessels compared to VEGF 164 tumors; 
however, whereas VEGF 164A108 _ 118 tumors are markedly 
more proliferative than VEGF 164 tumors [7], VEGF 188 
tumors show negligible growth [6]. We may conclude 
that this discrepancy indicates that the resulting vascular 
beds in the two tumors are different in ways that the 
intratumoral vessel density cannot measure, for example, 
in their connectivity to peritumoral vasculature [6], 
which may be a limiting factor for nutrient delivery; we 
propose that a separate aspect of VEGF signaling con- 
trols this latter behavior. Our model suggests a crucial 
difference between the two tumors: while VEGF 188 
tumors show minimal VEGFR2 activation, even less 
than that of VEGF 12 o tumors, VEGF 164A108 _ 118 tumors 
instead show even higher VEGFR2 activation than 
VEGF 164 tumors (Figure 8A). As a result, it may be pos- 
sible that the ability of the vasculature to support tissue 
growth may be dictated by VEGFR2 activation at the 
vessel itself, due to the receptor binding behavior of 
soluble VEGF, as suggested in numerous studies [16,17]. 
Furthermore, notice that peak VEGFR2 activation 
occurs between the VEGF 165 and VEGF 12 i isoforms 
(Figure 8Ai). This may support the observation that 
VEGF 12 o is at least as tumorigenic as VEGF 164 in some 
systems due to intrinsic differences between systems 
[25,39]. Overall, we hypothesize that intratumoral vessel 
density and branching is determined by matrix-bound 
VEGF detected by filopodia removed from the vessel 
surface, while vessel efficacy is dictated by VEGFR2 acti- 
vation at the vessel body. As a result, we thus posit that 
the tumorigenic behavior ultimately depends on the 
ability of VEGF to keep its receptor and NRP1 binding 
domains intact. The inability of VEGF 188 tumors to elicit 



VEGFR2 activation, despite it potentially having higher 
affinity to NRP1 than VEGF 164 , is a result of its degrada- 
tion and cleavage (Figure 7C). 

This hypothesis may also explain the paradoxical roles 
of specific proteases in tumorigenesis. Our mechanism of 
MMP-mediated VEGF redistribution (Figure 9A) shows 
that proteases will increase the functional soluble VEGF 
concentration by inhibiting the process of VEGF degra- 
dation. However, while several modes of VEGF redistri- 
bution, such as through MMP9, heparinases, or VEGF 
inhibitor cleavage are implicated as pro-tumorigenic, 
mediating an angiogenic switch [15,16,19,20,29,30,32], 
there are several notable exceptions where protease activ- 
ity instead leads to diminished tumor growth [7], an 
effect similar to the plasmin-mediated loss of wound 
healing due to a loss of angiogenesis [75]. We note an 
important trend: in studies where pro-tumorigenic beha- 
vior occurs from VEGF release, VEGF has not been 
shown to be directly cleaved. Bergers et al have shown 
that MMP9 mediates VEGF-release induced carcinogen- 
esis in pancreatic islets [16] without determining the 
mechanism of release. However, subsequent studies have 
shown that MMP9 does not necessarily cleave VEGF 
[15,31-33], as suggested by [7], but instead acts to cleave 
HSPGs directly [15]. In fact, Joyce et al [30] show that 
MMPs and heparinases have similar effects in the pan- 
creatic islet system, which strengthens the argument that 
the MMP9 induced angiogenic switch in [16] may have 
been mediated by HSPG cleavage. In contrast, studies 
where proteases reduce the angiogenic potential of VEGF 
show direct evidence of VEGF being cleaved and/or 
degraded [7,75,76]. We propose that VEGF needs to 
maintain coreceptor domains for effective tumorigenesis. 
Cleavage of VEGF, while increasing the total soluble 
VEGF concentration, may decrease overall VEGFR2 acti- 
vation due to a loss of NRP1 binding (recapitulated in 
Figure 8Aii); heparinases, MMP9, and VEGF inhibitor 
proteases also prevent VEGF degradation but redistribute 
intact, coreceptor-binding VEGF. 

Our results suggest that a central facet of VEGF pat- 
terning in vivo is its degradation, which we show neces- 
sarily occurs in an isoform-specific manner. Several 
mechanisms may underlie VEGF's isoform-specific 
degradation, and it is not currently known what, if any, 
mechanisms operate in the different in vivo experimen- 
tal systems used. Interstitial cells and endothelial cells 
from nearby vessels may uptake VEGF isoform in an 
HSPG- or NRP1- dependent manner. This is supported 
by observations that NRP1 and VEGF receptors are typi- 
cally present in many types of parenchyma, e.g. hind- 
brain [77], astrocytes [78], tumor [27,37], and skeletal 
muscle [79]. On the other hand, degradation may be 
due to the action of VEGF-degrading proteases, possibly 
the same proteases that also initially cleave VEGF; 
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however, this possibility does not seem to be consistent 
with developmental systems where VEGF cleavage has 
not been detected [17]. An interesting possibility that 
has recently been raised is that of VEGF inhibition by 
soluble VEGF inhibitors, e.g. sVEGFRl [11,33] or con- 
nective tissue growth factor [32], operating in an iso- 
form-dependent manner through either HSPG 
complexation or through NRP1 complexation. Interest- 
ingly, each of these mechanisms also supports the ability 
of proteases to allow VEGF to escape degradation, 
through either HSPG binding or NRP1 binding, mediat- 
ing VEGF redistribution. An important uncertainty in 
VEGF catabolism is whether endothelial cells represent 
the primary source of VEGF receptors and degradation 
in vivo or not. For example, while the background 
neural progenitor cells in the hindbrain may express 
NRP1 [78], endothelial expression is typically thought to 
be very strong [80]. Furthermore, immunochemical 
staining usually shows that VEGF concentrations dimin- 
ish precisely at the vascular front [12,13]. 

Besides being able to reproduce experimental observa- 
tions of VEGF patterning in vivo, is there evidence for 
isoform-specific degradation in tissues? VEGF (specifi- 
cally, VEGF bioactivity) has been shown to degrade in 
vitro under cell culture conditions [53,81] and in fact, 
proteases that cleave VEGF into shorter isoforms also 
seem to degrade it further [51,82]. In addition, while 
VEGF degradation has not typically been studied as a 
cause of VEGF patterning, this view is commonly 
accepted in numerous developmental systems (e.g. 
Decapentaplegic, Wingless in Drosophilia) [52,83,84]. 
Evidence for isoform-specific degradation however may 
come from an indirect source. Perlecan knockdown in 
zebrafish establishes a diffusible VEGF phenotype [85]. 
Surprisingly, it also increased total tissue VEGF levels. 
Relevant to the present study is the fact that total VEGF 
levels did not decrease, supporting the view that HSPGs 
may be important mediators of VEGF degradation. 
However, there are other pieces of evidence that seem 
to contradict isoform-specific degradation: intrave- 
nously-injected bevacizumab shows significantly greater 
tumoral deposition in VEGF 189 -expressing tumors com- 
pared to VEGF 165 tumors and VEGF 12 i tumors [86]. 
Preliminary computational results suggest that this 
observation supports the HSPG-binding-only model on 
the basis of the total number of bevacizumab binding 
sites, i.e. VEGF, increases proportionally to the VEGF 
isoform matrix affinity. 

An interesting prediction of the isoform-specific 
degradation model is that the extracellular residence 
times of different VEGF isoforms are roughly equal. 
Thus, a test of the model can be made by injecting 
labeled VEGF isoforms into tissues and measuring their 



half-lives in the interstitial fluid or lymph. If this test 
shows that the longer isoforms have significantly greater 
residence times in tissues, it would disprove the iso- 
form-specific degradation model. The residence time of 
VEGF in tissues is a direct measure of the overall degra- 
dation and clearance rate, and its constancy is identical 
to the statement that the total levels of VEGF in tissue 
are roughly constant against differences in patterns of 
isoform secretion and VEGF proteolysis, a finding that 
is suggested by results in [7,19]. Note that this state- 
ment, however, does not contradict isoform-specific 
degradation. In this model, the rate of degradation (or 
clearance) of the soluble fraction of VEGF is isoform 
specific, with heavier isoforms showing more rapid 
degradation; however, accounting for all phases of 
VEGF (e.g. matrix bound, receptor bound), the average 
rate of degradation of any VEGF isoform is nearly iden- 
tical (refer to Additional file 1, section S2). 

Overall, our results form the basis for a different view of 
VEGF patterning and endothelial behavior in response to 
VEGF. The assumption behind how VEGF patterning is 
intuitively interpreted is that of the transient: transiently, 
MMPs elicit VEGF release, which can increase VEGF 
receptor signaling on endothelial cells [16], and HSPGs do 
hinder diffusion, forming isoform-specific differences in 
soluble VEGF patterning. However, this assumption 
ignores what happens to the VEGF distribution over much 
longer periods of time, which are likely just as important 
for slowly evolving processes such as vascular patterning. 
The transient assumption only seems valid in studying in 
vitro systems, systems where VEGF is cleared very slowly 
[66]. On the other hand, in vivo systems seem to represent 
a major phenomenological difference due to their much 
more rapid VEGF dynamics (x < 1 h, Appendix Al), neces- 
sitating a steady-state analysis. Several additional assump- 
tions have also been made in our analysis (Table 1). For 
example, we assumed that the intrinsic proteolytic cleavage 
rate of all isoforms is identical. However, experimental stu- 
dies indicate that VEGF 189 may be more resistant to 
MMPs than VEGF 165 [7] . This point is interesting to note 
since another study, [51], found that VEGF ni (a form of 
VEGF also resistant to degradation or cleavage) also results 
in angiogenesis with high vascular density [51], similar to 
VEGF 164A108 _ 118 formed vessels [7]. Since the inability to be 
further cleaved/degraded may be the common theme, it 
may indicate that an inability to be cleaved/degraded alter- 
nately underlies higher vascular densities in those systems. 
Finally, in the current study, we specifically tested HSPG 
binding as the mechanism of isoform specificity in VEGF 
degradation. However, our results do not change if the iso- 
form-specific degradation occurs solely through soluble 
VEGF being degraded in an isoform-specific manner {not 
shown), such as in a NRPl-dependent fashion. 
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Conclusions 

In the present study, we have identified a general 
mechanism of VEGF transport that explains experimen- 
tal data regarding VEGF isoform patterning and proteo- 
lytic release at steady state, that of isoform specific 
degradation. We have further highlighted possible 
mechanisms by which information in the VEGF distri- 
bution can be used to guide vascular patterning in an 
attempt to explain vascular branching complexity and 
also the regulation of angiogenesis and tumor growth by 
proteases. A major limitation of our present study is 
that we were only able to broadly consider isoform 
ordering with respect to one or two features of the 
VEGF distribution. Instead, sprout formation and the 
subsequent patterning may be the result of a complex, 
temporal orchestration of multiple extracellular VEGF 
fractions, receptor signaling states, and cell types. Such 
an analysis is outside the scope of our study, however 
rule-based cell models and other modeling efforts are 
being developed to the sophistication required to 
address the multifactorial problem of vascular pattern- 
ing, e.g. [87-89]. Many studies implement VEGF degra- 
dation and/or ECM binding and proteolysis [90-93] 
however not in a form, that as we show, gives rise to 
differential VEGF gradients seen in vivo [13]. We believe 
that the advances made in our study, specifically the 
necessity of a mechanism of isoform-specific degrada- 
tion, would be useful to such studies to improve the 
predictions of the VEGF distribution. 

Important biological questions remain to be addressed, 
especially regarding the existence and nature of isoform- 
specific VEGF degradation, the nature of the MMP 
induced angiogenic switch, and the possible role of soluble 
inhibitors such as sVEGFRl in sprout formation, isoform 
patterning, and in the protease-mediated angiogenic 
switch. 

Additional material 



Additional file 1: Supplementary Methods and Results. The 

supplemental material is used to provide parameter estimation of VEGF/ 
HSPG binding affinities, provide in-depth methodology for the 
simulations, and derive important theoretical results. (File is in PDF 
format, readable with Adobe Reader). 
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